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Abstract — The present paper is focused on linear quadratic regulator (LQR) synthesis. It is evident that the solution of the 
problem leads to finding a solution of symmetric root locus (SRL) problem. It is proven and analytical evidences is also 
presented that by using SRL synthesis method and further application of LQR method guarantees stability, minimum phase 
and gain margins of the closed loop system. The presented mathematical relationships are further proven with experimental 
investigation in MATLAB programming environment. 
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I. Introduction 

The scope of the present paper is LQR synthesis in the state space. According to the state feedback control chosen, the 
synthesis procedure is related to determination of the state feedback matrix elements in a manner that a compromise between 
achieving the desired settling time and the control effort needed. In the presence of immeasurable states, observers are 
typically used for state vector estimation. The prerequisites needed for the synthesis are the system to be completely 
controllable and observable although advanced state variable design techniques can handle situations wherein the system is 
only stabilizable and detectable. In the present paper it is assumed that the system is completely controllable and observable 
and all the state variables are measurable. 

The procedure for LQR synthesis requires the feedback matrix to be chosen in order to satisfy the condition for minimum 
value of the cost function presented as integral performance estimation. The closed loop is considered optimal as it 
guarantees the minimum possible value of the cost function. The synthesis procedure allows the operator to choose the 
weights so that the compromise between the state cost and the control cost is acceptable. The major disadvantage of this 
technique is that the exact poles position and the control effort value are initially unknown. 

The present paper aims to investigate the relationship between the linear quadratic regulator procedure and the symmetric 
root locus of the closed loop system. 

Lor the provision of this analysis few main aspects should be investigated in details and further proven: 

1. Find analytic solution of the LQR synthesis based on the Pontriagin’s minimum principle. 

2. Show that LQR synthesis can be solved by SRL. 

3. Use SLR procedure to prove that when LQR synthesis is applied the closed loop system will always be stable and 
minimum phase. 

II. LQR SYNTHESIS BASED ON THE PONTRIAGIN’S MINIMUM PRINCIPLE 

Lor a given closed loop system, denoted in Lig.l, the state space open loop system model is: 

x = Ax + Bu 

y = cx (1) 
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The feedback control law in the presence of a reference signal is: 

u(r) = -Kx(r)+Gv(r) ^ 

where v(f) is the reference value and the matrix G is given an auxiliary function to recalculate the reference value for steady 

state error free system [1], Assuming the following transformation v(?) = 0 t > 0 the control design problem is known as 

regulator problem [2], It has the following aspects: determination of K so that all initial conditions are derived to zero using 
a specified approach determined by the desired system performance which corresponds to the closed loop pole values. 

The optimal control design (Fig.l) is a constrained optimization problem with cost function J (x,u, t, t f ) which can be 
expressed as: 


J = 


(f/ )Sx(f/ ) + i [ 7 , *Q X + 7 , u ' Ru ]^ 


(3) 


and the constrain given by equation (1) [4,6,7]. The matrices in (3) must satisfy S > 0, Q > 0 and R > 0. The implementation 
of co-state vector /. [5] consisted of the Lagrange multipliers transforms (3) into unconstrained optimization problem, 
described by: 


/ = ^x r (t / ) s x(t / ) + J 


^x 7 Qx + ^u r Ru + X T (Ax + Bu - x) 


dt 


(4) 


In order to obtain the minimum value of (4) all the derivatives of J with respect to x, u and x(t f j must be equal to zero. 
According to the last operation the closed loop dynamics can be determined by the following simultaneous equations for a 
given mixed boundary conditions x(0) = x 0 and ~t-[t f ) = Sx (t f j [4]: 


x = Ax + Bu 
i = - Qx-a t i 

u=RB^ (5 ) 

This is the so called two point boundary value problem (TPBVP) [8]. The system (5) can be presented as a homogenous 
vector-matrix differential equation: 


x 

i 


A 

-Q 


BR B r 

-A 7 


Where H is the so called Hamiltonian matrix [5], 
The substitution: 

X = Px 


( 6 ) 


(7) 
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and the equation system (6) leads to the following expression: 

i. = Px + Px = Px + PAx - PER B r Px = -Qx - A r Px 

After the reduction of the common multiplier x equation (8) can be written in the form: 

P = -A r P - PA - Q + PER B r P 

known as differential vector-matrix Riccati equation. 

This equation can be solved backward in time under the condition: 

p( v )=s 

It defines the optimal control law: 

u(r) = -R'B r P(r)x(r) = -K(r)x(r) 

The steady state value of equation (11) is the solution of the algebraic Riccati equation: 

A r P + PA + Q - PER B r P = 0 
and determines the control law in the form: 

u(f) = -R ‘B r Px(r) = -Kx(f) 

III. LQR SYNTHESIS AND SRL 


( 8 ) 

(9) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


The closed loop system dynamics is described by the equation (6). The eigenvalues of matrix H are the closed loop poles. 
They are doubled in number than the number of the eigenvalues of matrix A and can be further determined as a solution of 
the following equation: 

det (si - H) = 0 (14) 

The implication of the substitution Q = C C and introduction of the following notations are further used for the purpose of 
this analysis: 


W P (s) = C(s/ - A) ‘ B 
W,', (-s) = B r (-H-A r )"‘c 7 
/3(s) = det(.yI-A) 


(15) 

(16) 
(17) 

According to some well known matrix relations and transformations described in [7] the equation given below can be 
obtained: 


det(jI-H) = (-l)"£)(j)£)(-j)det[l + R- 1 W;(-j)Wp(j)] 

where n is the state vector dimension. 

For SISO system equation (18) can be presented as: 

det(vI-H) = (-1)" £>(^)/)(-j)[l + R- 1 Wp (s)W p (-.v)] = 0 
The open loop system transfer function could be then described as a polynomial ratio: 
. . . . _i N(s) 

W P (s) = C(sl —A) B = -^ 


(18) 


(19) 


( 20 ) 
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Considering (20) equation (19) is transformed into: 

det(,I-H) = (-l)"[/3(v)/?(-,)H-^ 1 fV(v)fV(- J )] = 0 

Consequently the closed loop poles determination can be solved as standard root locus problem with variation of R ' into the 
range (0,oo) . The specifics here are that the locus is symmetric in accordance with the abscissa as well as with the ordinate 

and has 2 n loci. Half of them are positioned in the left half of the complex plane and are related to the n stable closed loop 
poles corresponding to the state vector x. The other half of loci are in the right half plane and are related to the n unstable 
closed loop poles which correspond to the co-state vector k. The unstable poles can be transformed into stable if the 
differential Riccati equation is solved backward in time[7,8]. 

IV. Stable and minimum-phase closed loop system 


As it was described above the weight matrices must be Q > 0 and R > 0. It is also known that only the state cost weight can 
be equal to zero. Thus according to the root locus it is more convenient weight factor p to be implemented under the 
assumption of R = 1 [2]. The cost function can be then described by the following integral relation: 

1 r oo 1 fco 

J = — px Qx + u u dt <=> = — py y + u u dt 

2 Jo (22) 

The SRL is determined by the equation [6]: 

D{s)D(-s) + pN(s)N(-s) = 0 (23) 

and is expresses as a function of the weight factor p which varies into the range [0,co). 

For p = 0 minimum settling time is not a requirement. The stable closed loop poles coincide with the system roots in the left 
half plane and/or with the reflection of the unstable open loop poles in accordance with ordinate. Thus the closed loop system 
stability is guaranteed with lower values of the control signal. 

For higher values of p the poles are shifted following the root loci and respectively the control effort increases too. 


For /? —»co, m closed loop poles approaching open loop zeros in the left half plane and/or their reflection in the right half 
plane. In this so described case the impact of those poles and zeros is negligible. The rest (n-ni) poles approaching infinity 

following specific trajectories and form the so called Butterworth pattern. The settling time is shortest possible however the 
control effort rises dramatically which in some cases is unacceptable [2,8]. 

The solution of the following equation is investigated in order to confirm those statements: 

A W = D (s)D(-s) + pN ( s)N(-s ) = 0 (24) 

with variations of p within the range [0, go) . In accordance with the below polynomial form: 


D^s) = s" + d 1 s n ~ 1 H-1- d n 

and 

N(s) = n 0 s m + « l s m_1 4-h n m = n 0 N i (s) 

is obvious that: 

N(s)N(-s) = (s)^ (-s) 

For p = 0 : 


(25) 


(26) 


(27) 


A(s) = D(s)D(-j) = 0 (2g) 

Thus the roots of A (A), which are closed loop poles, coincide with the open loop poles and their ordinate symmetric 
reflections. 
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For p —> oo the first addend in equation (24) l)(s) D(-s) is neglected: 

A(i) = pN^s)N(—s) = pnlN { (s)N l (—s) = 0 


(29) 


The polynomial A(s) is sum of the two addends D(s) D(~s) and pipN l (s ) A) (- 5 ) . Since the zeros of /V, (.v) /V, (-,v) result 

in zero sum and zero second addend as well then they respectively zero the first addend too. Consequently it can be presented 
as: 


(30) 


D(s)D(-s) = N\ (s)^ {~s)D t ( s)D l (s) 
and A(s-) is further transformed into the equation: 

A(-y) = N, (s)N t {-s)D t ( s)D t ( -s) + pn 2 0 N 1 ( s)N , (- 5 ) = 

= (s)D 1 (- 5 ) + pnl ] = 0 

A iW M») (31) 

The first multiplier Aj ( 5 ) roots are the zeros of the open loop system and their symmetric ordinate reflections. The second 
multiplier: 

A 2 ( J ) = D l ( A )' D 1 {~ S ) + P n 0 
can be transformed as follows: 


(32) 


At ( v) = (—1)" s 2( " + d t s 


2(n-,n)-l 


+ d 2 s 


2(n-m)-2 


' + ^2 (n-m) + P n 0 


For p —> 00 the roots of A, (.v) converge to the roots of equation: 

(-!)—/—)+ p n i = o 

Respectively the roots of A, (.sj converge to the roots of the binomial equation: 


(33) 


(34) 


p {n-m) 


- pn„ =0 

\n—m ' u 


(' ^ (35) 

The roots of equation (35) appear to be the apexes of regular (2(n-m)J - dimensional polygon and are positioned on the 

/ 2 \l/2(n-m) 

circle centered in the origin with radius I pip 1 .For higher values of p they are shifted from the origin following the 

trajectories by the circle radiuses forming equal angles ——— . The resulted plot with (2 (n-m)) rays is associated with 

\n — m) v ; 

the so called Butterworth pattern [1,2,3], 

V. Numerical example 

The open loop system investigated is unstable and non minimum-phase described by the below presented transfer function 

[4]: 


W(s) = - 


(s-2)(s-4) 


s 2 (s-1)(s-3)(s 2 +0.8s + 4.16) 

The problem formulated is LQR synthesis aiming to minimize the value of the cost function / (22). 

Specific algorithm has been developed in MATLAB for the purpose of this analysis. The resulted SRL is shown in Fig.2. 


It is evident that: 
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> tor p = 0 two of the SRL loci begin at the two stable complex roots of the open loop system s l 2 = -4 ± 2 i . Another 
pair loci begin at zero pair roots ,S ' 34 = 0 and third pair begin at the ordinate symmetric unstable roots ,v. = I and 
s 6 = 3 . The analysis shows that in this case in particular the closed loop system is at least considered stable; 

> with the increase of the value of p the last pair loci described above are approaching the ordinate symmetric positive 
zeros s = 2 and s = 4 ;when p approaches infinity that pair ends at those zeros which leads to compensation of the 
influence of the zeros s = — 2 and s = -4 with two of the poles. The rest four loci are shifting from the origin 
following the radiuses of central circle, forming equal angles of nj4. The result is retention of the closed loop system 
stability, short settling time, and increase in the natural system frequency; however the system requires higher control 
signal values. 


Root Locus 



VI. Conclusion 

The paper presents analytical evidence that LQR synthesis can be considered as SRL problem. 

In the presence of evidence for completely controllable and observable (stabilizable and detectable) open loop system the 
state feedback synthesis is invariant in accordance with open loop system stability. 

LQR can be also taken into consideration as a method for state feedback matrix generation with guarantees for stability, 
minimum phase and stability margins of the closed loop system. 

This type of analysis, using graphic interpretation (ex. Fig.2) allows direct identification of the damping ratio and 
corresponding natural frequencies co n and respectively further calculation of the weight coefficient p = cop for any point 
positioned on the loci. 

The results delivered give further option for investigation of possible relationship and application of the pole placement 
synthesis and LQR synthesis based on SRL in the frequency domain. 
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